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ABSTRACT 

Three-dimensional Doppler tomography has been used to study the Ha emission sources in the RS 
Vulpcculac interacting binary. The 2D tomogram of this binary suggested that most of the emission 
arose from the cool mass losing star with additional evidence of gas flowing close to the predicted 
trajectory. However, the 3D tomogram revealed surprising evidence of a more pronounced gas stream 
flow at high V z velocities from -240 to -360 km s _1 . This behavior is most likely caused by magnetic 
activity on the cool star since the central velocity plane, defined by V z = km s" 1 , should be coincident 
with the orbital plane of the binary if the flow is dominated by gravitational forces only. RS Vul has 
been detected as both an X-ray and a radio source, and it is possible that the RS Vul gas stream may 
have been deflected by magnetic field lines. This flow is distinctly different from that found in the 
streamlike state of U CrB, in which the gas flow was confined mostly to the central velocity plane. 
Subject headings: techniques: image processing - accretion, accretion disks - stars: binaries: close - 

(stars:) binaries: eclipsing - (stars:) circumstellar matter - stars: imaging - stars: 

individual (RS Vulpeculae,U Coronae Borealis) 



1. INTRODUCTION 



The extension from 2D to 3D tomography is now a 
reality. Two-dimensional Doppler images can be cre- 
ated with the standard Filtered Back Projection (FBP) 



Method {e.g., Kaitchuck et al. 1994) by projecting the 



sources of emission onto a single central (V x ,V y ) velocity 
plane. This plane is expected to coincide with the orbital 
plane of the binary since this is the plane in which the 
gravitational forces dominate the gas flows. However, 
these sources may have non-zero velocities beyond the 
central 2D plane since they are not geometrically thin 
and are not confined to the orbital plane. Instead, 3D 
images can be derived over a range of V z values transverse 
to the central plane. The resulting 3D maps are consis- 
tent with the 2D Doppler tomograms obtained with the 
FBP method and reveal substantial flow structures be- 
yond the central velocity plane. The 3D reconstruction 
method applied in this work is the Radioastron o mical 
Approach (RA) developed by Agafonov ( |2004a|b[ ) ; Aga- 
fonov & Sharova ( 2005a|b I . A detailed discussion of the 
differences betw een this technique and the FBP method 
can be found in |Agafonov, R ichards & Sharova ( 2006 1 
and Agafonov, Sharova fe Richards | ( I2009p . In particu- 
lar, the RA method is especially useful when the number 
of projections is relatively small because it can achieve 
the same resolu tion as FBP while using fewer profile s 



(see Figure 1 of Agafonov, Richards & Sharova 



2006) 



The effectiveness of the RA method has been examined 
by using initial 3D test models, calculating ID profiles 
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from t hese models, an d then performing the reconstruc- 
tion. ( jSharovaT 2006 ) compared the computed models 
with the initial models and showed that the computed 
emission maxima coincided with the initial model, and 
all details were reproduced. So, there was very good 
agreement between the initial test model and the recon- 
structed 3D tomogram. 

The two-dim ensional Doppler tomograms of the A lgol- 
type binaries ( Richards, Albright & Bowles 1995 Al- 
bright fc Richards ||1996| |Richards ||2004l display a, aT 



verse range of circumstellar structures. Ihese include a 
gas stream, accretion annulus, transient accretion disk, 
shock regions and sometimes a chromospheric emission 
source in the short-period (P < 5 days) Algols, as well 
as a classical ac cretion disk in the long-perio d (P > 5 - 
6 days) Algols. [Richards fc Albright I ( ji999| ) found that 
some systems, like U Coronae Borealis (U CrB) and U 
Sagittae (U Sge), alternate between streamlike and disk- 
like states and display changes that can occur within a 
year, and sometimes overnight. They identified RS Vul 
as a member of this group. 

The first applications of the 3D procedure focused on 
the interacting Algol- type binary U CrB. The 2D tomo- 
grams constructed from Ha spectra of this binary col- 
lected in 1993 and 1994 displayed two distinct patterns: 
a transient accretion disk representing the disklike state 
(in epoch 1993) and a bright gas stream flowing along the 
predicted ballistic trajectory representing the streamlike 
state (in epoch 1994). The 3D tomograms of these emis- 
sion sources were even more distinct than the 2D ver- 
sions, and both states showed evidence of a prominent 
emission source associated with the mass gainer, result- 
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ing from the impa ct of the gas stream onto the photo 
sphere of that star flAgafonov, Richards fc Sharova |2006 



|Agafonov, Sharova fe Ricliarcls""|PooTf| r Chromospheric 
emission from the mass loser was also visible especially 
when the binary was in the disklike state. The 3D im- 
age also showed distinct evidence of gas flows (jets) with 
high V z velocities. These jets were associated with three 
emission sources: (1) the impact/splash region where the 
gas stream strikes the mass gainer in the streamlike state 
(V z ~ -300 to +100 km s _1 ); (2) the interaction between 
the gas stream and accretion disk in the disklike state 
(V z ~ 300 to 400 km s" 1 ); and (3) the Localized Re- 
gion where the disk material interacts with the incoming 
gas stream after traveling around the mass gainer (V z 
h200 to +500 km s" 1 ). These high values arc ex- 
pected since the star-stream impact occurs at speeds of 
up ^500-600 km s ~\ compared to the slower rotation 
of the mass gainer (Richards 1992; Blondin^ Richards fc 



Malinkowski ]1995| [Richards fc Ratliff |1998[ ). Moreover, 
the pattern of V? velocities associated with the disk emis- 



patte 

sion implies that the disk may be inclined to the orbital 
plane; which subsequently suggests that the rotation axis 
of the mass gainer could be precessing as a consequence 
of the star-stream impact. This result would not have 
been possible from 2D tomograms alone. 

In this paper, we have continued our exploration of 
3D tomography by studying the RS Vul system. Specifi- 
cally, we have 1) identified the features in the 3D images 
and related them to those seen in the 2D tomograms; 2) 
compared the 3D tomogram of RS Vul with those of U 
CrB; and 3) provided a physical interpretation of these 
features. The system parameters and 2D results are de- 
scribed in Section 2, the 3D tomograms are described in 
Section 3, a model for RS Vul is given in Section 4, and 
the Conclusions are outlined in Section 5. 

2. 2D DOPPLER TOMOGRAPHY OF RS VUL 

The Algol-type binary RS Vulpcculac is an interact- 
ing system in which the cooler Gl III secondary star 
is transferring mass to its B5 main-sequence companion 
(the primary) through Roche lobe overflow. RS Vul is 
a shor t-period Algol (P or b — 4.4776635 da ys; |Holmgren| 
11 1989 1 wit h an orbital inclination of 78.7° ( |Hutchings &;| 
Hi ll |1971 ) . The properties of the system are M p = 6.59 



Mq, M s = 1.76 M Q , R p = 4.71 R Q , R s = 5.84 R G , sys- 
temic velocity V — -20.1 km s _1 , velocity semi- amplitude 
K v = 54.0 km s _1 , with a mass ratio, q=0.27 (Holmgren 
19891. Eighty-one Ha (6562.8 A) spectra with a reso- 



lution of 0.166 A/pixel or 7.6 km s 1 were collected in 
1993 with the KPNO 0.9m Coude Feed Telescope (see 



Richards fc Albright | [T999). The spectra were collected 



at closely spaced positions around the entire orbit of the 
binary. Since the observed spectra are dominated by 
the spectrum of the B5 mass gaining star, model atmo- 
spheres calculations were used to represent the photo- 
spheric contributions of the stars, and these stellar con- 
tributions were subtracted fr om the observed spectra to 



create difference spectra (see Richards fc Albright [1999 



for details). The difference spectra were used to calculate 
the tomograms. 

The standard 2D Doppler tomogram of RS Vul 1993 is 
shown in Figure 1 along wit h the 2D tomogra m of U CrB 
in the streamlike state (from |Richards ( 2001[ )). These 2D 
images were compared with the 3D Doppler tomograms 



based on the same data. In the tomograms, the solid 
trajectory is the gravitational free-fall path of the gas 
stream; and the circles along this trajectory are marked 
at intervals of a tenth of the distance from the Li point to 
the distance of closest approach to the mass gainer. The 
largest solid circle and the smaller dashed circle mark the 
inner and outer edge of a Keplerian disk, respectively; the 
asterisk is the predicted location where the gas stream 
should strike the photosphere of the mass gainer; and the 
plus sign marks the center of mass of the binary. The 
RS Vul tomogram shows the presence of the gas stream, 
a circumprimary emission source, chromospheric emis- 
sion, and a regio n dominated by absorption called the 
absorption zone (Richards 20011. The strongest emis- 
sion source corresponds to the c hromospheric Ha emis- 
sion ( |Richards fc Albright |[l996| . 

3. 3D DOPPLER TOMOGRAPHY 

3.1. Reconstruction of the 3D Tomogram 

For the reconstruction of the 3D Doppler tomograms, 
we anal yzed the same se t of Ha fine profiles that were 
used by Richards (2001) to construct the standard 2D 
Doppler tomogram of RS Vul. The 3D Doppler tomo- 
gram was calculated as a set of numerical values in the 
cells of a cube, with dimensions (V x ,V y ,V z ) having val- 
ues ranging from -700 to +700 km s . The maximum 
of this function was normalized to the unit to permit 
the visualization of the image based on any slice taken 
from any direction. The velocity resolution of the image 
depends on the number of projections (i.e., number of 
spectra) and their distribution in orbital phase; and the 
orbital inclination influences the ratio of the resolutions 
in the V x , V y , and V z directions. The wavelength resolu- 
tion of the spectra is 7.6 km s _1 , so it plays a smaller role 
in setting the resolution of the image since it is already 
sufficient to produce the best possible velocity resolution 
for the tomogram. In particular, the velocity resolutions 
in the (V x ,V y ), (V Z ,V X ), or (V z ,V y ) planes are equal only 
for an inclination of 45°; and the resolutions in the V x 
and V y directions are better than that in the V z direction 
for high values of inclination. Since the orbital inclina- 
tion of RS Vul is high (i — 78.7°), the 3D tomogram 
was restored with a resolution of 30 km s _1 in the V x 
and V y directions and 105 km s _1 in V z direction. Since 
the resolution in the V z direction is approximately 3.5 
times lower than in the other directions, features on the 
tomogram may look elongated in this direction. The ap- 
pearance of the images has been simplified by displaying 
the (V x ,V y ) slices in the horizontal plane to be consis- 
tent with the direction of the orbital plane of the binary, 
and the slices containing the V z axis are in the vertical 
direction. 

A direct comparison was made between the observed 
spectra and those computed from the reconstructed 3D 
Doppler tomograms. Figure 2 shows this match for the 
3D Doppler tomograms of RS Vul. The radial velocity, 
V r is plotted versus orbital phase, with a resolution of 30 
km s _1 . The observed spectra are displayed in the left 
frames and the spectra computed from the 3D Doppler 
tomograms are displayed in the middle frames. Two S- 
wave patterns are noticeable in the observed and com- 
puted spectra corresponding to the motions of the two 
stars. The overall agreement between the observed and 
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computed spectra is very good. The quality of the fit is 
shown in the right frame of Figure 2, which displays the 
relative chi-square statistic, x 2 /xli versus orbital phase. 
Here, x 2 is normalized to the critical value of x 2 j which 
corresponds to the largest acceptable value of x 2 at the 
99% confidence level (e.g., |Skilling fc Bryan |1 984). Fig- 
ure 2 shows that x 2 was less than 20% of the critical value 
at most orbital phases and it was only as high as 65% of 
the critical value over a narrow phase range from 0.47 to 
0.49, perhaps because the 3D code assumes that the gas 
is optically thin when it might be optically thick at these 
phases. So, in all cases, x 2 was much lower than the crit- 
ical value which confirms that the 3D reconstruction has 
produced a reliable computed model for RS Vul. 

3.2. Description and Interpretation of Emission 
Features 

RS Vul and U CrB have similar geometries, so we ex- 
pect to identify similar features in their Doppler images. 
Figure 3 displays 20 slices in the (V x ,V y ) plane for val- 
ues of V z ranging from —540 km s _1 to +540 km s _1 , 
in 60 km s _1 steps. The superposition of the color or 
gray-scale images with the contour images allows us to 
emphasize the main features in each V z slice of the 3D 
image. Figure 4 displays the images of seven representa- 
tive V z slices over the range from V z = —480 km s _1 to 
+480 km s _1 . A similar figure for U CrB in the stream- 
like state is shown in Figure 5 for comparison. Figure 6 
displays several cross-sections of the flow in the (V y ,V z ) 
plane. In Figures 3-6, the light (yellow) areas with blue 
contours represent absorption while the dark (deep pink) 
a reas with red contours represent emissi on . 



km s^ 1 in Figures 3 and 6 (see also Table 1). Moreover, 
it is distinctly different from the same feature seen in U 
CrB 1993, which displayed a symmetric V z distrib ution 
from V z = -180 to +180 km s" 1 (see Table 1 in |Aga- 
fonov, Sharova & Richards 2009[ ). The circumprimary 



Agafonov, Sharova & Richards (2009) identified sev- 
eral features in the 3D tomogram of U CrB that were 
consistent with t hose found in the 2D image studied by 
Richards (20011: (1) circumprimary emission (centered 
on the mass gamer) and an accretion annulus (a circular 
region of low velocity) ; (2) chromospheric emission from 
the cool donor star (the secondary); (3) the gas stream 
flowing along its predicted ballistic trajectory; (4) the 
star-stream impact region, where the gas stream strikes 
the stellar surface; (5) emission within the predicted locus 
of the accretion disk; (6) a localized region between the 
stars where the gas stream strikes material that has cir- 
cled the mass gainer; (7) a second localized region where 
the gas stream makes impact with the outer edge of the 
accretion disk; (8) a jet or flare associated with the chro- 
mosphere of the secondary star; and (9) a high-velocity 
stream or jet extending in the Vz-direction beyond the 
stream-star impact region. 

The 3D image of RS Vul displays all of the emission 
sources found in U CrB (see Figs. 3 and 4), except that 
there was almost no emission from an accretion disk and 
little evidence of the star-stream impact region. Figure 3 
shows that the strongest emission sources are associated 
with the cool mass loser as well as the mass gainer. The 
emission sources are described in Table 1 along with their 
typical velocity ranges and a brief description of each 
feature is provided below. 

1. The circumprimary emission centered on the velocity 
of the mass gainer in RS Vul is labeled (1) in Figure 4. 
It is the strongest source of emission in the 3D image 
(see Figure 3, slice at V z = —160 km s _1 ). This feature 
can be seen in the slices with V z velocities from to -320 



emission has higher velocities than expected from syn- 
chronous rotation, so the star may have been spun up by 
the impact of the gas from the donor star. 

2. Chromospheric emission from the cool mass losing 
star is feature (2) and it is comparable in strength to 
the circumprimary emission. It is associated with a wide 
range of V z velocities from —530 to +350 km s _1 (see 
Figures 3, 4, 6, and Table 1). This range is illustrated 
in Figure 6 and resembles a jet-like feature in the V z 
direction. The feature has the same velocity as the mass 
losing star. 

3. The gas stream flowing along its predicted ballistic 
trajectory is feature (3). It flows from the LI point to- 
wards the mass gaining star with speeds up to 450 km s _1 
in the (V x ,V y ) plane, with a V z velocity of —300 km s _1 . 
RS Vul is unusual because the flow along the gas stream 
is not strongest in the V z = km s _1 frame but in the 
V z = -240 to -360 km s" 1 range (see Figs. 3 and 4). This 
is distinctly different from U CrB in its streamlike state 
(Fig. 5), in which the gas stream flow was strongest in 
the central velocity plane (V z = km s _1 ) and located 
directly on the predicted trajectory of the gas stream in 
the (V x ,V y ) plane. This is not the case for RS Vul. 

4. The predicted location of the star-stream impact re- 
gion, where the gas stream strikes the stellar surface, is 
found at (4) near the location where the gas stream tra- 
jectory intersects the inner part of the disk. Even in its 
most extended state (when it is at V z — -300 km s^ 1 ), 
the gas stream (Feature 3) never gets to this predicted 
velocity, but is truncated at the velocity corresponding to 
that of Keplcrian flow at the surface of the mass gainer. 

5. There is very little evidence of any emission within the 
predicted locus of the accretion disk (Feature 5) in the RS 
Vul images, although some clumping is seen within this 
region. This feature was much stronger in U CrB, and 
there is slightly more space around the mass gainer in 
U CrB to form a disk. However, gas may have traveled 
around the mass gainer and created the localized regions 
(described below) when the interaction with the stream 
slowed the gas flow. 

6. A Localized region (LR) between the stars labeled (6) 
was found near the center of the tomogram where the 
flow around the mass gainer is expected to strike the in- 
ner ed ge of incoming gas strea m (see Car tesian models of 
.gPerj Richards ||1992[ [l993| , RW Tau ( |Vesper fc Hem" 1 
eycutt |l993p , and other direct-impact systems). In RS 



Vul, this region is detected at high positive V z velocities 
ranging from +50 to +300 km s _1 and seems to develop 
into the structure seen at the LI point from V z = +360 
to 540 km s _1 (Fig. 3). The classical localized region 
between the stars is seen along the V x —0 km s _1 line for 
V z velocities from 50 to 540 km s _1 in Figure 3. How- 
ever, it is difficult to say whether the part with LI point 
velocities in the (V x ,V y ) plane for V z = +180 to +300 
km s _1 (Fig. 4) should be considered part of the cool 
star or part of the LR. This feature is also seen in Figure 
6 slices for V x — —20 to +20 km s _1 . This kind of braking 
was first suggested by Richards (1992) in the study of the 
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Algol-type binary (3 Persei and confirmed using hydrody- 
namic simulations of Blondin, Richards & Malinkowski 
H Q1995[ ) and IRichards fc Ratliff | (l998p . Similar results 
were found later fr om g as dynamical rnodel i ng of cata- 
clys mic variables by Kuznet sov et al. 1(|2001[ ), |Bisikalo et 



aL| ( |2000a| ), and jfjisikalo et al I ( |2000b7 
T. A second compact region labeled (7) also has high 
V z velocities from to 350 km s _1 and is located on the 
V y = km s _1 line (see Figures 3 and 4). In Figure 6, 
it can be seen in the V x — 100 to 120 km s _1 frames, 
with an average V z velocity of about 250 km s _1 . This 
feature may later merge with the other localized region 
that runs along the line of centers between the stars. 

8. The feature labeled (8) has a velocity similar to the 
V x and V y velocities of the donor star, but with high 
positive V z velocities in the range from Vj=+180 to +540 
km s^ 1 (see Figures 3 and 4). It has V x and V y velocities 
of 100 to 300 km s" 1 (see Figs. 3, 4, and 6). This 
feature could be the result of chromospheric emission on 
the cool mass loser, and the high velocities suggest that 
it could correspond to a jet or flare. However, another 
explanation is that it is related to the interaction between 
the inner edge of the gas stream and the material that 
has circled the mass gainer. A similar emission source 
was found in U CrB. 

9. The absorption feature in RS Vul seen in the lower 
left and right quadrants of Figure 3 from V z = -120 to 
+ 540 km s -1 was called the absorption zone by Richards 

[] pool " 



2004), and is associated with the locus where 



the gas temperature becomes too high to emit at Ha. 
A structure at a similar location was identified in the 
ultrav iolet tomogram of U Sge by |Kempner fc Richards 



(119991 and in the Ha tomogram of U CrB by Agafonov, 
|SIiarova fc Richards | ( [2009| ). 

The 3D tomograms also reveal an interesting symmetry 
in the emission sources found in the V z slices. A strong 
source near the mass losing star (in velocity space) is de- 
tected with intensities of 0.852 and 0.955 in symmetrical 
locations at V z — —300 and +300 km s _1 , respectively 
(see Figure 4). A similar symmetric intensity pattern 
(I = 0.343 and 0.333) is seen at V z = -480 and +480 
km s _1 . In the case of the mass gaining star, another 
symmetric pattern is revealed. The most intense emis- 
sion feature in the V z — 180 km s _1 frame in Figure 4, 
with 1=0.612, is located along the V x = km s _1 line 
at V y — +80 km s _1 , in the opposite V y direction from 
the circumprimary emission. These symmetries may be 
purely coincidental or perhaps we can see the same fea- 
ture with motions on both sides relative to the central ve- 
locity plane, with a —V z component in — z direction and 
with a +V Z component in +z direction. Nevertheless, it 
is expected that the gas flows should be symmetric about 
the central velocity plane if gravitational forces dominate 
the gas flows. So, the displacements of emission features 
(1), (2), and (3) from the central velocity plane should 
not be unusual. 

4. A MODEL FOR RS VUL 

Figure 7 displays the 3D tomograms of RS Vul and U 
CrB in a 3D velocity cube together with the 2D inte- 
gral images for the main planes: (V x ,V y ), (V X ,V Z ), and 
(V y ,V z ). These 2D integral images represent the inte- 
grals along the V z , V y , and V x directions, respectively; 
and the image in the (V Xl V y ) plane is equivalent to the 



2D Doppler tomogram. The U CrB images in this figure 
are derived from the combined disk and stream states of 
that system. The 3D velocity cube of RS Vul illustrates 
that RS Vul has no gas stream emission along the pre- 
dicted trajectory in its central velocity plane, although 
the 2D integral image in the (V x ,V y ) plane shows that 
the gas stream emission is indeed prominent, but at V z 
velocities much higher than expected. In addition, both 
binaries display a gas velocity distribution that is inclined 
to the vertical direction (V x = km s _1 ) in the 2D inte- 
gral (V X ,V Z ) plane, but vertical in the 2D integral (V y ,V z ) 
plane. 

A pattern has emerged from the study of the 3D to- 
mogram of RS Vul. Both the circumprimary emission 
source labeled earlier as (1) and the gas stream (3) were 
offset toward negative V z velocities, while the chromo- 
spheric emission (2) had V z velocities over the range from 
—530 to +350 km s _1 (see also Figure 6). So, the cir- 
cumprimary emission was strong over the same range of 
V z velocities corresponding to the strongest gas stream 
emission. Moreover, the gas stream emission extended 
to the Keplerian velocity at the photosphere of the mass 
gainer, while the mass gainer was spinning faster than 
the synchronous value. Therefore, the results are consis- 
tent with the assumption that the circumprimary emis- 
sion (from a photospheric bulge) is produced when the 
gas stream strikes the surface of the star. In addition, 
the mostly negative V z velocities for the gas stream and 
the circumprimary emission suggest that the gas does not 
move in the central velocity plane, as expected if grav- 
itational forces dominate the gas flows, and so the flow 
around the mass gainer is concentrated on one side of the 
central velocity plane. 

While the image reconstruction process can now pro- 
duce 3D velocity images, we are still unable to create the 
desired 2D or 3D Cartesian images because the various 
emission features have different velocity fields. However, 
an estimate of the gas stream flow can be made from 
information provided in Figure 1. Based on the point 
where the gas stream is truncated, we can produce a 
Cartesian model to show the dominant emission sources 
in RS Vul. The distance from the Li point to the center 
of the primary star has been divided into 10 equal seg- 
ments in the center frame of Figure 1, and the similar 
small circles on the Doppler tomograms are marked at 
velocity intervals corresponding to those distance inter- 
vals. Comparing the center and right frames of Figure 1, 
we see that the gas stream flow is truncated after trav- 
eling 7 small circles on the gas stream trajectory, which 
is equivalent to 70% of the distance from the Li point 
to the center of the primary star. The star-stream in- 
teraction occurs along the path of the gas stream right 
where the gas stream ends (at 70% of the distance from 
the Li point). Both Figures 3 and 4 show that the gas 
flow reaches the stellar surface corresponding to the solid 
circle in the tomogram, which corresponds to the Keple- 
rian velocity at the surface of the mass gainer. Beyond 
this point, there is no evidence of any significant disk in 
this binary at this epoch. 

Evidence of circumstellar gas distributed well above 
and below the orbital plane of the binary has been de- 
rived from studies of ultraviolet resonance lines and the 
Ha line of several Algol systems (e.g., Peters and Pol- 
idan 1984, Richards 1993). This gas can reach heights 
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comparable to the radius of the mass gaining star. How- 
ever, semi-analytical ballistic studies (e.g., Lubow and 
Shu 1975, 1976) or hy drodynamic simulations of the mass 



transfer process (e.g., Richards & Ratliff 19981 assume 



that the gas should be found very close to the orbital 
plane, if mag netic fields are not involved. However, [Stern] 



ct al. 



( 1992 1 suggested that coronal mass ejections could 
account for more than 10% of the mass transfer expected 
solely from Roche lobe overflow by gravitational pro- 
cesses alone. The coronal behavior of the cool stars in Al- 
gol binaries is enhanced because these stars rotate more 
rapidly than the slowly rotating Sun. If we can scale 
the solar coronal behavior to the enhanced levels seen 
in many Algol binaries, then radio and X-ray flares de- 
tected from these systems can be inte rpreted as the peri - 
odic reconnection of the coronal field flStern et al. |1992 1 
RS Vul was d etected as an X-ray source by | White fe Mar 



shall 



(1983| with an X-ray luminosity, L x = 2.0 x 10 JU 
erg s _1 , which is about half of the x-ray luminosity of 
(3 Per. Radio emission from RS Vul was also detected 



by Umana et al (1998) with a flux density of 0.26 mJy 



and radio luminosity, L ra dio — 3.7 x 10 16 ergs cm 
Hz -1 , compared to 0.30 mJy and 9.1 x 10 16 ergs cm -2 
s _1 Hz -1 for U CrB. So, the magnetic fields of the active 
stars in RS Vul and U CrB should be comparable. 

The long-term radio flare survey of the p rototype Algol 
system (3 Per by Richards et al. | ( |2003| ) demonstrated 
that major flares (up to Uy) occur regularly and pre- 
dictably every 48.9 ±1.7 days (or 17 orbital cycles), thus 
the flaring rate is at least 7.5 yr _1 and could be as high as 
100 yr _1 (for weaker flares every orbit). Mass ejections 
from X-ray flares observed with EXOSAT and GINGA 
are es timated to be ~ 5 - 50 xl0~ 15 M 
1992), and a total mass ejection rate of : 
x 10 4 Mq yr _1 would correspond to the detected flar- 
ing rate of 7.5 yr _1 for j3 Per. These arguments suggest 
that coronal mass ejections may influence the transfer of 



(Stern et al. 
3.75 - 37.5 



gas between stars and even the angle at which the flow is 
directed. The differences between the gas distributions 
in RS Vul and U CrB cannot be readily explained by 
differences in magnetic activity based on their x-ray and 
radio luminosities. However, it is evident that magnetic 
fields may play a role in the distribution of circumstcllar 
gas in both systems. 

5. CONCLUSIONS 

Three-dimensional Doppler tomography has been used 
to study the Ha emission sources in the RS Vulpeculae 
interacting binary and has led once again to the discovery 
of gas flows with significant V z velocities. The 2D tomo- 
gram of this binary suggested that most of the emission 
arose from the cool mass losing star with additional ev- 
idence of gas flowing close to the predicted trajectory. 
However, the 3D tomogram revealed surprising evidence 
that a more pronounced flow of gas along the stream was 
not found as expected in the central velocity plane (V z 
— km s -1 ) but at V z velocities from -240 to -360 km 
s _1 . Magnetic activity on the cool mass losing star may 
explain the deflection of the gas stream flow away from 
the orbital plane of the binary where the gravitational 
forces dominate the interactions between the stars. This 
flow behavior is distinctly different from that found in the 
streamlike state of U CrB, in which the gas flow was con- 
fined mostly to the central velocity plane. Exploration of 
the 3D velocity structures of other binary systems should 
improve our understanding of these out-of-plane gas mo- 
tions. 
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Fig. 1. — Two-dimensional Doppler tomogram of RS Vul - 1993 (left) taken from |Richards | ( |2001[ ) , a Cartesian model of RS Vul (center), 
and the 2D tomogram of U CrB in the streamlike state (right). In the tomograms, the solid trajectory is the gravitational free-fall path of 
the gas stream; and the circles along this trajectory are marked at intervals of a tenth of the distance from the Li point to the distance of 
closest approach to the mass gainer. The largest solid circle and the smaller dashed circle mark the inner and outer edge of a Keplerian 
disk, respectively; the asterisk is the predicted location where the gas stream should strike the photosphere of the mass gainer; and the 
plus sign marks the center of mass of the binary. 
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Fig. 2. — Comparison between the original data (left frame) and the spectra computed from the reconstructed 3D Doppler map of RS 
Vul (middle frame) in terms of the radial velocity, V r versus orbital phase. The right frame displays the orbital phase variation of the 
relative chi-square statistic, x 2 /Xcj where \c ' s * ne critical value corresponding to the 99% confidence level. The agreement between the 
observed and computed spectra is very good. 
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Fig. 3. — (V x ,V y ) slices of the 3D Doppler tomogram of RS Vul, for velocities beyond the central velocity plane from V z = -540 km s 1 
to V z = +540 km s _1 , in steps of 60 km s _1 . The images are shown for intensity levels: +0.01, 0.05, 0.1, 0.3, 0.5, 0.7, 0.9 (red) and -0.01, 
-0.05, -0.1, -0.3 (blue). The central source, corresponding to the mass gainer, is strongest at V z = - 160 km s — 1 and very weak at Vz= 
km s , while the gas stream emission is strongest from V z = -360 km s — 1 to V z = -240 km s — 1 . 
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Fig. 4. — Visualization of the seven most interesting (V x ,V y ) slices in the 3D Dopplcr Tomogram of RS Vul (1993) displayed symmetrically. 
The main features are: (1) circumprimary emission, (2) chromosphcric emission from the donor star, (3) gas stream, (4) star-stream impact 
site (the asterisk *), (5) the predicted locus of accretion disk, (6) localized region between stars, (7) other localized region, and (8) high 
intensity region with velocity close to that of the donor star. 
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Fig. 5. — Visualization of the seven most interesting (V x ,V y ) slices in the 3D Doppler Tomogram of U CrB (1994) in the streamlike 
state displayed symmetrically. The main features are: (1) circumprimary emission, (2) chromos p heric emission from the donor star, (3 ) 
gas stream, (10) a high-velocity jet in the Vz direction (sec Agafonov, Richards & Sharova (2006); Agafonov, Sharova & Richards (2009). 
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Fig. 6.— Cross sections of the RS Vul 3D Doppler tomogram in the (V y ,V z ) plane for V x = -380, -200, -120, -80, -20, 0, 20, 100, 120, 200 
km s — 1 with intensity levels of 0.01; 0.05; 0.10; 0.30; 0.50; 0.70; 0.90. The numbered features correspond to those in Figure 3. The maps 
for V x values from +200 to +380 km s _1 contain no information. 



Three-dimensional Doppler Images of RS Vul 



11 




RS Vul (P = 4.48 d ) 



U CrB (P = 3.45 d ) 



E o 

5 




Vx (km s" 1 ) 



it i < 





■200 200 400 

Vx (tons" 1 ) 




-600 - 400 -200 200 400 

Vy (km s ') 



Vx (km S ) 



2D integral images 

Vx,Vy 
(integral along Vz direction) 



Vx,Vz 
(integral along Vy direction) 

■^Central velocity plane^* 




-000 -400 -200 



Vy,Vz 
(integral along Vx direction) 

*-Central velocity plane-* f 




Fig. 7. — The 3D tomograms of RS Vul and U CrB displayed in a 3D velocity cube (top row) together with the 2D integral images for 
the (V x ,V y ) plane (second row), (V X ,V Z ) plane (third row), and (V y ,V z ) plane (last row). The 2D integral image in the (V x ,V y ) plane is 
equivalent to the 2D Doppler tomogram. 
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TABLE 1 

Characteristics and locations of prominent emission features in the 3D tomogram 







Location: central velocity 


or velocity range (km s ) 


No. 


Emission Feature 


V x 


Vy 


Vz 


1 

2 
3 
4 
5 
6 
7 
8 


Circumprimary emission 
Chromosphcric emission 
Gas stream flow 
Star-stream impact region 
Locus of the accretion disk 

Localized Region (LR) - Part 1 (between stars) 
Localized Region — Part 2 
Jet or flare near donor star 


(-60 to +60) 
(-120 to +120) 

-520 to a 
no source found 
very weak 

100 

200 (100 to 300) 


-60 (0 to -120) 
210 (60 to 360) 
100 to 300 a 

to 100 
-80 to 30 
~200 


-160 (0 to -320) 
-160 (-530 to +350) 
-300 (-500 to -60) 

180 (50 to 300), 450 (360 to 540) 
240 (150 to 350) 
360 (180 to 540) 



a corresponds to the predicted ballistic trajectory of the gas stream. 



